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ABSTRACT 



Studies of beryllium abundance in large samples of solar-type stars show a small fraction of extremely beryllium- 
deficient stars, which challenges our current understanding of light element depletion in these stars. We suggest two 
possible scenarios that may explain this high level of Be depletion: early accretion and rotational mixing. We show that 
in both cases, the conditions required to reach the observed level of Be depletion are quite extreme, which explains 
the very small fraction of detected Be outliers. We suggest that substantial Be depletion can be obtained in stars if 
they were fast rotators in the past, with high initial rotational velocities and short disc lifetimes. Our analysis suggests 
that rotational mixing may not be efficient enough to deplete Be in less than 10 Myr. Consequently, the detection of 
strongly Be-deficient stars in clusters younger than ~ 10 Myr may provide a genuine signature of accretion process 
and the proof that some protostars may undergo many extreme bursts of accretion during their embedded phases of 
evolution. 
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^ 1. Introduction 

> ■ 

^ 8 The light chemical elements lithium, beryllium and boron 
, . provide good tracers of internal mixing processes in stel- 
lar interiors. A huge amount of observational and theoret- 
ical works has been devoted to the analysis of their abun- 
dances at the surface of various types of stars. Interest 
in these elements has recently increased within the con- 
text of exoplanet discoveries because of possible links be- 
tween Li depletion in a star and the presence of a planet 
orbiting this star (see e.g Sousa et al. 2010 and refer- 
ences therein). This suggestion motivated extensive stud- 
ies of solar-type stars with the main goal of analysing 
Li and Be in planet hosting stars and of confirming 
or refuting the above-m e ntioned link dTa keda et al] 1201 it 
iGalvez-Ortiz et a.1,1 120TTL iDekado Mena et al.ll2oTa . This 
has yielded large samples of homogeneous data, which pro- 
vide excellent opportunities to study properties and trends 
of light-element depletion. One striking property resulting 
from these studies is a small fraction of highly Be-deficient 
stars in sufficiently large samples, which questions our cur- 
rent understanding of Be destruction in solar- type and low- 
mass stars. 

In a sample of 82 stars. [Santos et aTl (|2004f ) found three 
abnormally Be-depleted objects, namely HD 4391, HD 
20766 and HD 280807. These Be-deficient stars are dwarfs 
with T e ff in the range 5700 - 5900 K. Lithium depletion 
is poorly constrained but is also significant, with Li abun- 
dances A(Li)Q < 1 (upper limit). According to the measured 
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vsini of these stars, their rotation rate is similar to that 
of th e Sun, and th e ir e stimated ages are between 1 and 2 
Gyr. iTakeda et al.l (|201ll ) found four strongly Be-depleted 
objects within their sample of 118 stars: HIP 17336, HIP 
32673, HIP 64150, and HIP 756 76. These obje c ts hav e sim- 
ilar properties as those of the ISantos et al.l ((2004) sam- 
ple: strong Li depletion, an age of ~ 3 — 5 Gyr, and low 
vsini ~ 2 km/s. The properties of these seven objects are 
summarised in Table [U Existing calculations for Be deple- 
tion hi_thJi5j^ on rotational mixing scenar- 
ios ([Pi nsonneault et al. 1989) or mixing by internal waves 
(Montalban & Schatzman 2000 ), cannot explain these low 
Be abundances (see e.g Fig. 4 oflDelgado Mena et al.ll2012l ). 



1 The notation A(E) defines the abundance of element E and 
is equal to \og[N (E) / N (H)] + 12, where N(E) and N(H) are the 
mole fraction of element E and of hydrogen, respectively. 



In this paper, we explore two scenarios that may ex- 
plain these Be-deficient outliers, which are observed in 
an effective temperature range 5000 — 6000 K and corre- 
spond to a mass range 0.8 — 1.1 M . The first scenario 
is based on episodic accretion at early stages of evolution 
(©? which can strongly affect the internal structure of 
young accreting objects and enh ance the deple t ion o f light 
elements, following the id eas of iBaraffe et al.l (|2009l ) and 
iBaraffe fc Chabrierl (|2010f ). The second scenario is based 
on rotational mixing. We adopt a simplified treatment of 
rotational mixing in stellar interiors, which we describe and 
validate in ^3l We quantitatively show (Q that substantial 
Be depletion can be obtained in stars if they were fast rota- 
tors in the past. This idea was mentioned in Tak eda et al.l 
([201lh . but no quantitative studies exist to confirm it. 
Discussion and conclusion follow in Sect. [5j 
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Table 1. Parameters of strongly Be-depleted objects. , , , , , 



Star 


^eff 


A(Be) 


A(Li) 


Age 


v sin i 




(K) 






(Gyr) 


(km/s) 


HD 4391 (a) 


5878 


0.64 


< 1.09 


1.2 


2.72 


HD 20766 (a) 


5733 


< -0.09 


< 0.97 


1.6 


1.98 


HD 20807 (a) 


5843 


0.36 


< 1.07 


2.1 


1.74 


HIP 17336 (b) 


5671 


< -0.85 


< 0.8 


3.89 


2.07 


HIP 32673 (b) 


5724 


< -0.78 


< 1.0 


3.16 


2.83 


HIP 64150 (b) 


5800 


< -0.88 


< 1.0 


4.2 


2.21 


HIP 75676 (6) 


5772 


< -0.96 


< 1.1 


4.16 


2.18 



Note s. (a ) From iSantos et al.1 (j2004f ) (b) From iTakeda et all 
(pOlOh and ITakeda et al.1 (|2011l ). 



2. Episodic accretion scenario 

iBaraffe et all (|2009l ) and iBaraffe fc Chabrierl (|2010l ) re- 
cently investigated the effect of episodic accretion at very 
early stages of evolutio n on the structure a nd Li deple- 
tion of low-mass objects. IBaraffe et al.l (|2009l ) showed that 
bursts of accretion can produce objects significantly more 
compact than their non-accreting counterpart of same mass 
and age. These authors suggested that non-steady accre- 
tion during the embedded phase of protostar evolution pro- 
duces the observed spread in luminosity in the Herzsprung- 
Russel diagram (HRD) of youn g clusters. As a consequence 
of the more compact structure, IBaraffe fc Chabrierl (|2010f ) 
showed that objects undergoing very strong accretion 
bursts have significantly higher central temperatures and 
can be severely Li-depleted. This scenario provides an ex- 
planation for the unexpected level of Li depletion observed 
in a few low-mass members of young clu sters (~ a few Myr). 
Based on the same calculations as in Baraff e~fc Chabrierl 
(|2010l) . we find that significant Be depletion can take place 
in stars that have undergone short and very intense bursts 
of accretion with M ~ 5 10 -4 M yr _1 . Those results are 
illustrated in Fig. [TJ which displays the abundance of Be 
as a function of time at early stages of evolution of low- 
/solar-mass stars. These calculations show that the early 
history of accretion can significantly deplete Be in stars 
in the mass range of interest. Note that Li is totally de- 
stroyed in the mode l s dis played in Fig. [TJ As explained in 
IBaraffe fc Chabrierl (|2010l ), the more compact and hotter 
structure of accreting models increases the maximum tem- 
perature reached at the bottom of the convective envelope, 
which increases the level of Li and Be destruction com- 
pared to the non-accreting counterparts of same mass. We 
find that substantial depletion of Be only happens if the 
protostar undergoes several bursts (more than ten) of ex- 
treme intensities with rates ^> 10 -4 M yr _1 . According 
to numerical simulations of collapsing cloud cores, typical 
mass accretion rates onto the pr otostar vary in the range 
10~ 6 M^yr" 1 - 1Q ~ 4 M^r" 1 (jVorobvov fc Basul 120101 : 
iMachida et al.l 12010). Several burst events with rates ex- 
ceeding a few times 10 -4 M yr _1 up to ~ 10 -3 M yr _1 are 
predicted, but only for high clou d masses with high initia l 
angular momentum values (see IVorobvov fc Basul l2010l ). 
This suggests that high levels of Be depletion by accre- 
tion process should only be produced under quite extreme 
conditions and should therefore be rare events. 
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Fig. 1. Abundance of Be (divided by the initial abun- 
dance Beo) versus time (in yr) in stars that have under- 
gone a given number A/burst of accretion bursts with M=5 
10" 4 M yr _1 and duration At burst = 100 yr. The blue 
solid line corresponds to an accreting sequence with final 
mass 1 Mq with 7Vb U rst=20 bursts. The red dashed line 
corresponds to a final mass 0.8 M with iVb U rst=16. The 
accretion process starts from an initial seed mass of 10 Mj 
(see Baraffe & Chabrier 2010 for details). The black dotted 
line indicates the Be abundance for non-accreting models 
of 0.8 Mq or 1 M , which are the same for both masses 
since Be is not depleted in those models. 



3. Rotational mixing scenario: the model 

In the following, we explore quantitatively whether rota- 
tion can explain extreme Be depletion, based on the im- 
plementation of rotation effects in the same stellar evo- 
lution code as in the previous section to study accretion 
effects (see Baraffe et al. 2009 and references therein). 
Different parametrisations for the treatment of rotational 
effects in ID st ellar evolution code s exis ts in the lit- 
erature, see e.g. IPinsonneault et all (ll989t): iHeger et all 
(2000); iDenissenkov fc Toutl (|2QQQh : iMaeder fc Mevnetl 
(2010) and references therein. These formalisms are usu- 
ally based on the t reatment of lEndal fc Sofia! (Il978l) (e.g . 
IPinsonneault et al.l [1989) or on the mode l of IZahnl (11 992) 
and its subsequent developments (see e.g. IMaeder fc Zahnl 
1998). All these formalisms contain free parameters that 
are calibrated to reproduce observations (e.g solar surface 
rotational velocity, surface abundances of light elements) 
and yield very different results when applied to the sam e 
case study, as recently illustrated by iPotter et al.l (|2012f ) . 
Strong uncertainties still remain regarding the treatment 
of rotation and its effect on stellar structure and evolution, 
despite very sophisticated and complicated formalisms. As 
our aim here is to explore a possible connection between ro- 
tation and extreme depletion of light elements, we chose to 
follow a simple approach to explore the effects of rotation 
on Li/Be depletion. 
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3.1. Rotation law 

Our approach adopts the hypothesis of iBouvier et al.l 
(|l997l ). namely: 

1. Solid body rotation - Instead of solving an angular mo- 
mentum equation, we assume solid-body rotation. The 
underlying hypothesis is that angular momentum is ef- 
ficiently transported within the interior. The physical 
processes transporting angular momentum can be of 
magnetic origin (e.g. magnetic torque) and/or hydrody- 
namic origin (shear turbulence, large-scale currents, in- 
ternal waves). The efficiencies and relative contributions 
of these processes to the total angular momentum trans- 
port are still poorly understood. Note that helioseismol- 
ogy shows that the bulk of the sol ar radiative zone is i n 
nearly solid-body rotation, see e.g. lGarcfa et al.l (|2007l ). 

2. Disc locking - In the early stages of evolution, the star is 
efficiently coupled to its accretion disc and the rotation 
rate Q remains constant during the disc lifetime: tt(t < 
Td) fto. 

3. Magnetic breaking - Angular momentum is lost through 
magnetic breaking by the stellar wind. The rate of loss is 
computed using Kawaler's law (|Kawalerl ll988). includ- 
ing saturation of the losses when Q > Q s &t- We use the 
parametrisation 

o / /? \ 1/2 / 7M" \ -1/2 

'm b = -tfn»(-) (— ) ifn^ (i) 

and 

J mh = -Kn 2 sat n(—) (_) otherwise. (2) 

We use ^ sat =14 ^© (|Bouvier et al.lll997l ). K is a con- 
stant that is calibrated to yield the correct surface ve- 
locity v ~ 2 km/s of a solar model at t = 4.6 Gyr. 
This yields a value K = 2.7 x 1Q 47 , in agre ement with 
previous works (see e.g. IBouvier et al.lll997l ). Magnetic 
breaking is applied from the end of the disc locking 
phase. 

Within this approach, the rotational history of a star 
is entirely determined by two free parameters: and Td- 
Values adopted for these parameters are discussed in Sect. 

3.2. Rotational mixing of chemical species 

Assuming that transport of angular momentum is efficient 
enough to maintain solid-body rotation, we only need to 
treat rotational mixing of chemical species. Rotational mix- 
ing in stably stratified regions driven by dynamical flows 
(e.g. meridional circulation, shear instabilities) is imple- 
mented in stellar evolution models to explain various ob- 
servational features (see e.g Maeder & Meynet 2010 and 
references therein). It is taken into account by introduc- 
ing an effective "turbulent" diffusivity D t in the evolution 
equation of chemical species, 

dXi id, 2 ox v 
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Fig. 2. Evolution of the rotation rate for an initially slow 
(lower red curves) and fast (upper blue curves) rotating 
solar-mass star, for different disc lifetimes: Td = 1 Myr 
(solid curves); Td — 5 Myr (dashed curves); Td = 10 Myr 
(dotted curves). The solar value is indicated by a solid cir- 
cle. Direct and inver ted triangles are observational data 
compiled in Table 1 of lBouvierl (|2008f ) and represent respec- 
tively the 10th and 75th percentiles of observed rotational 
period distribution. 

The expression for D t significantly varies from one au- 
thor to the other (see references above and also Potter et 
al. 2012). We chose to adopt a simple prescription (see e.g. 
iRichard et al.ll2004l ) that will permit us to capture the es- 
sential effects of rotational mixing: 

D t = or\U(r)l (4) 

where a is a free parameter and U (r) is related to the radial 
velocity component of the meridional circulation: 

u r (r,0) = U(r)P 2 (cos0) J (5) 

where P 2 is the second Legendre polynomial. In the solid- 
body rotation hypothesis, U (r) is given by (Richard et al. 
2004) 

U(r) = 1 En, (6) 

pgTc p M* V ad - V 

where Eq is defined by 
and M* is given by 

with p m the mean density within the sphere of radius r. 

We neglect the effect of molecular weight gradients, 
which have stabilising effects and consequently decrease the 
efficiency of mixing. Molecular weight gradients appear in 
stellar interiors when a significant amount of central hydro- 
gen has been burned during the main sequence evolution. 
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Fig. 3. Calibration of a for different rotational histories (represented with the same colour-coding as in Fig. [2]) of a 1 
M star using solar values for the ratio of present to initial abundances of Li (Li/Lio=l/140, see left panel). The black 
error bar shows the range 100 < d(Li) < 200. The right panel shows the corresponding depletion of beryllium and the 
calibrated values of a. In both panels, the case without rotational mixing is shown for reference (black continuous line). 



Since most of the rotationally induced depletion occurs ear- 
lier, i.e. for t < 0.1 Gyr, before any significant /i-gradient 
appears in the interior (see ^3.3|) . neglecting /i-gradient ef- 
fects in our analysis will not change the conclusions of this 
work. 

This model is simple to implement in a stellar evolution 
code and has the advantage of involving only one free pa- 
rameter a that needs to be calibrated (see ^3.3|) . We stress 
that in this model, the rotation history is totally decoupled 
from the mixing process, i.e. a has no feedback on Q(t). 

3.3. Rotational histories and calibration of a 

We adopted diffe rent rotational histories based on the work 
of lBouvie"rl (|2008[ ) . We used three values for the disc lifetime 
Td'. 1 Myr, 5 Myr and 10 Myr. These values are represen- 
tative of disc fraction analysis in young clusters of different 
ages, with a disc fr action exponent ially d ecaying betwee n 
- 1 and - 10 Myr (|Mamaiekll2009h . As in lBouvier] (|2QQ8h . 
we constructed rotational models for slow and fast rotators, 
with initial rotation periods Po equal to 10 days (2.5 f£ ) 
and 1.2 days (~ 21 fi©), respectively. 

The different values of and Po mentioned above de- 
fine six different rotational histories. Figure [2] shows the 
evolution of the rotation rat e for each case . Observational 
data compiled in Table 1 of iBouvie r (2008) are shown for 
comparison. The models cannot reproduce the lower en- 
velope of observed rotational periods, failing at predicting 
sufficiently low rotation periods at t ~ 40 Myr (see Fig. 
[2]). This may stem from the solid-body rotation assump- 
tion. iBouvierl ([2008) was able to reproduce the lower en- 
velope of the distribution by assuming strong decoupling 
between the core and the envelope. An envelope rotating 
ten times slower than the core yield s the ob s erved low sur- 
face velocities. On the other hand, IBouvie 2 (|2008l ) showed 
that strong coupling between core and envelope, resulting in 
nearly solid-body rotation, reproduces the upper envelope 
of the velocity distribution. Therefore, solid-body rotation 
assumption is a good proxy for fast rotating models during 



the main-sequence (it is not excluded, however, that signif- 
icant differential rotatio n can develop during the PMS, see 
lEggenberger et al.1l2012f ). 

The calibration of the free parameter a defining the dif- 
fusion coefficient D t (see Eq. [4j) is based on the observed 
solar depletion of li thium, taken to be Li/Lio = 1/140 (see 
iRichard et al .11 19961 ). A 1 M star model is evolved with the 
six different rotational histories defined above and for each 
cases, a is calibrated to obtain the observed Li depletion 
at the solar age (£© = 4.6 Gyr). The result of this pro- 
cedure is shown in Fig. 03 It shows that beryllium is also 
depleted, with d(Be) ~ —0.2 dex at the solar age. Note 
that recent measurements of abundances in the solar pho- 
tos phere do not suggest a significant depletion in beryllium 
fsee lGrevesse fc Sauvallll998l : lAspfund et alJ2QQ9h . The cal- 
ibration process yields values of a in a range between 0.18 
and 5.30, as summarised in Fig. 03 Higher values of a are 
needed for lower rotation rates to reproduce the expected 
lithium depletion. The left panel of Fig. [3] shows that the 
depletion history of lithium is quite different for fast ro- 
tators and slow rotators: fast rotators have essentially de- 
pleted their lithium after t ~ 100 Myr and the depletion 
on the main-sequence is modest due to the low values of 
a, whereas slow rotators undergo a more steady depletion 
due to higher values of a. As expected, in our approach de- 
pletion increases with increasing a and increasing rotation 
rates. 

3.4. Validation of the model for different stellar masses/ages 

To additionally test the validity of our model assumptions, 
we explored the behaviour of the calibration obtained in 
£13.31 for different masses in the range 0.8 — 1.1M . Each 
stellar mass was evolved with the six adopted rotational 
histories (see Fig. [2]) and their corresponding value of a de- 
fined in £ 13.31 (see right panel of Fig. [3]). Results are displayed 
in Fig. do at three different ages: t=1.7 Gyr, 4.6 Gyr and 

2 Initial abundances of Li and Be were set to the meteoritic 
values A(Li)=3.3 and A(Be) = 1.42 fsee Grevesse fe SauvaJl998h 
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Fig. 4. Lithium (left) and beryllium (right) depletion versus effective temperature for four different masses M = 0.8 M 
(red symbols), 0.9 M (green symbols), 1 M (blue symbols), 1.1 M (yellow symbols) at three different ages t=1.7 
Gyr (crosses), 4.6 Gyr (dots), 6 Gyr (squares). For each mass, at a given age, the results of the six models corresponding 
to different rotational histories are displayed (see text). Initial abundances are shown as horizonta l long-dashed lines, 
lithium abundances below -1 are plotted at -1. Th e four dashed lines corres pond to models A of iPinsonneault et al.l 
(1990). The dotted curve corresponds to model C3 of Pinsonneault et al. (1990), which predicts the highest Be depletion. 



6 Gyr. The models predict a decreasing abundance of Li 
with T e ff; below 5500 K, all models have A(Li) < 1, which 
is roughly the detection limit. This is i n qualitative agree- 
ment with observati ons, see e.g. Fig. 5 oflGalvez-Ortiz et al.l 
(120111 ) and Fig. 5 of Taked aTetall (120101) . For a given mass, 
the spread of symbols at t — 1.7 Gyr over ~ 2 dex reflects 
the spread seen in the left panel of Fig. [3l 

Beryllium abundances show a much lower dispersion, 
because Be depletion is much less sensitive to rotational 
mixing and age. Our models show a slight increase of A (Be) 
with decreasing T e ff for the mass range considered here. The 
reason is that for decreasing stellar mass, the convective en- 
velope reaches deeper layers and rotational mixing becomes 
less and less efficient as a result of the decrease in r\U(r)\ 
and hence in the mixing diffusivity D t (see Eq. [3]). 

Different authors disagre e on the observed trend of 
A(Be) with T eff . Figure 9 of iBoesgaard fc Krugler Hollekl 
(2009) suggests a significant increase of beryllium abun- 
dance with d ecreasing T e g in the range T e ff ~ 5500 — 6500 
K; Fig. 6 of iDekado Mena et al.l (|2012l ) (see also Fig. 4 
of iGalvez-Ortiz et al. l201lh exhibits a clear trend of de- 
creasing beryllium abundance with decreasi ng T e ff in the 
range Teff = 5500 - 6200 K; whereas Fig. 9 of iTakeda et~aD 
(|201ll ) exhibits a constant beryllium abundance in the 
range T e ff ~ 5500 — 6000 K. Given these opposite results, we 
cannot conclude on the validity of our approach to predict 
the observed trend of A(Be) with T e ff. 

We show in Fig. [4] models of IPinsonneault et al.l (|l990f ) 
for the depletion of Li and Be at an age of 1.7 Gyr, 
used b y [Galvez- Ortiz et al.l (|201ll ) and lDelgado Mena et al.l 
(|2Q12l ) for comparison with their observations. These mod- 
els agree very well with our calibration of the depletion of 
lithium, whereas the right panel of Fig. [4] shows that our 
parametrisation yields a lower Be depletion than that of 
IPinsonneault et al.l (|l990l ). 
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Fig. 5. Beryllium depletion at t — 1.7 Gyr for our fast- 
rotating model (Po = 1.2 d) with = 1 Myr (continuous 
red lines) and = 5 Myr (dashed blue lines), for all four 
masses M = 0.8, 0.9, 1,1.1 M . The symbols correspond to 
different values of a: 1 (dot), 2 (diamond), 3 (left triangle), 
4 (right triangle), 5 (down triangle). Open symbols are the 
objects of Table [TJ triangles stand for upper limits. The 
initial Be abundance is shown as a horizontal dashed line. 



4. Rotational mixing scenario: results 

The simplified treatment of rotational mixing adopted 
in this work has been validated in the previous section. 
Calibrating the free parameter a based on the Sun's prop- 
erties and on a given rotational history yields values for a 
in the range ~ 0.2 — 5.3. This calibration provides realistic 
predictions for Li and, to a lesser extent, Be depletion as 
a function of mass (or T e ff) in the mass range of interest 
0.8 — 1.1 Mq. The model can therefore be used to explore 
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Fig. 6. Abundances of Li (left panel) and Be (right panel), divided by initial abundance, versus time for the fast rotator 
model with ra = 1 Myr and Po = 1.2 d, with a=2. Results are displayed for M = 0.8 M (dashed line) and M = 1.1 M 
(solid line). 



conditions required for extreme Be depletion, adopting the 
range of calibrated a values but for other rotational histo- 
ries, to explain the objects listed in Table [U These have 
common properties: low values of vsinz, roughly at the 
present day solar rotation rate (2 km/s), and an estimated 
age older than 1 Gyr. As shown in Fig. [2j different rota- 
tional histories yield models that converge towards the same 
track after ~ 1 Gyr. Consequently, observed Be-deficient 
stars could have been faster rotators than the Sun in the 
past, even though they now have low vs'mi. Note that we 
ignore the rotational history of the Sun and whether it was 
a slow or fast rotator in the past. 

Given this unknown, we assumed two rotational histo- 
ries that produce the fastest rotators at early ages, (jd = 1 
Myr, P = 1-2 d) and (r d = 5 Myr, P = 1-2 d), and adopted 
higher values of a, but within the calibration range. The re- 
sulting Be depletion at t = 1.7 Gyr is shown in Fig. [5j 
The objects of Table [T] are also shown for comparison. The 
figure shows that for the fastest rotating model (r^ = 1 
Myr, Po = 1.2 d), a value of a = 2 yields strong deple- 
tion (A(Be)< -1) for 1 M and 1.1 M stars, i.e. in the 
T e ff range of the observed outliers. For a slightly slower ro- 
tating model (rd = 5 Myr, Po = 1.2 d), higher values of 
a are required to reach a similar level of Be depletion, but 
still within the calibration range. Slower rotational histories 
were not explored since they would require values of a out- 
side the calibration range. Interestingly enough, significant 
beryllium depletion (A(Be) ~ -1 or log Be/Beo ~ —2) is 
difficult to achieve for M < 1 M with our model, except 
for unreasonable values of a outside the calibration range. 

Variations of Li and Be abundance versus time in the 
fastest rotating model (r^ = 1 Myr, Po = 1.2 d) with a=2 
are displayed in Fig.[6]for the two extreme cases M=0.8 M 
and M=l.l M . For these fast rotators, Li is rapidly de- 
pleted, decreasing by a factor 5 after 5 Myr and a factor 
100 after 10 Myr, whereas Be depletion starts later in those 
models, at ages > 10 Myr. Even with the highest value of 
a within the calibration range (a = 5), Be is depleted by 
only a factor two within the first 10 Myr. 

3 Results for ages > 1 Gyr are similar. 



5. Conclusion 

We have presented two possible scenarios that can explain 
extremely Be-deficient solar-type stars. The first scenario is 
based on the effect of accretion at very early stages of evo- 
lution of the protostar. For very intense bursts, we showed 
that Be can be very rapidly destroyed, producing strongly 
depleted objects at ages of only a few Myr, which could ex- 
plain the Be-deficient solar- type stars listed in Table [U High 
levels of Be depletion, as observed, can only be reached if 
the protostar undergoes several bursts (more than ten) of 
extreme intensities (M ^> 10 -4 M yr _1 ). Those events are 
predicted by numerical s imulations of cloud collapse (see e.g 
IVorobvov fc Basull2010f h but only for massive and initially 
rapidly rotating clouds. This scenario should therefore pro- 
duce only very few Be-deficient stars, which agrees with the 
few observed numbers (typically 3-4 objects in a sample of 
- 80-100 stars, see ®- 

The other possible scenario is based on rotational mix- 
ing. We have explored those effects based on a simpli- 
fied treatment of rotation in a ID stellar evolution code. 
Rotational mixing was described in our approach by a sin- 
gle free parameter. Our model assumed solid-body rota- 
tion throughout the stellar evolution of both slow and rapid 
rotators. Slow rotators are though t to develop a high de- 
gree of differential rotation (see iBouvierl [2008) and it is 
not excluded that rapid rotators can develop significant 
differentia l rotation during the pre-main sequence phase 
(see lEggenberger et al.ll2012[ ). Differential rotation induces 
shear mixing, which is by definition absent from our model. 
This is embedded in our single free parameter a, calibrated 
to reproduce the solar lithium depletion. We showed quan- 
titatively that a model assuming solid-body rotation that 
also has values of the free parameter within the calibration 
range can reach the observed level of Be-depletion assuming 
initially fast rotators with short disc lifetimes (r^ ~ 1 Myr). 
The observed level of Be depletion in the above-mentioned 
outliers requires initially very fast rotators with very short 
coupling times between the star and its disc. From a statis- 
tical point of view, those objects should belong to the upper 
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envelope of the observed rotational period distribution as a 
function of age and should th us only represe nt a small frac- 
tion of objects (< 10%, see e.g lBouvie r 2008). This could be 
consistent with the small fraction of observed Be-deficient 
solar- type stars. This exploratory study highlights the pos- 
sible connection between rotation and extreme depletion of 
light elements. Given our assumptions, we cannot rule out 
the possibility for slow rotators, which are thought to de- 
velop a high degree of differential rotation, to reach extreme 
depletion through an efficient mixing by shear turbulence. 
Additionnal studies should be devoted to investigate this 
possibility and confirm our findings. Ideally, these should be 
based on improved modelling of rotational mixing, includ- 
ing a consistent picture of angular momentum transport 
and turbulent mixing in stellar interiors. Improvement is ex- 
pected from the development of multi-D numerical simula- 
tions of stellar interiors, based either on existing numerica l 
tools ([Espinosa Lara fc Rieutordll2007t [Ballot e t al.l 120071 ). 
or on new tools (see e.g I Viallet et aT I I201lh . Multi-D nu- 
merical simulations will likely yield a better calibration of 
the free parameters inherent to current rotation formalisms 
used in ID stellar evolution codes and provide better un- 
derstanding of the different mechanisms leading to angular 
momentum transport and chemical mixing in stars. 

If accretion and rotation are two possible processes that 
can produce Be-deficient stars, it may be difficult to find 
different observational signatures that would allow distin- 
guishing one process from the other. Regarding a possible 
signature of accretion process, our analysis suggests that 
rotation is hardly able to deplete Be within the first ~ 
10 Myr of evolution of solar-type stars. It may be argued 
that differential rotation could enhance the depletion in this 
early phas e through additional m ixing by shear instabili- 
ties (see Egg enberger et all l2Q12h . In our 1.1 M model, 
the radiative core develops at £ ~ 1.4 Myr, and Be-burning 
starts at the centre of the star at t ~ 3 — 4 Myr. Given 
the size of the radiative core, we estimate that the mix- 
ing diffusivity should be enhanced by more than one order 
of magnitude, reaching values of about 10 7 cm 2 /s, to al- 
ter the abundance in the convective envelope, and thus at 
the surface, in less than t = 10 Myr. Although this de- 
pends on how strong differential rotation can develop dur- 
ing this early phase, it seems unlikely to reach such high 
diffusivities. Detection of strongly Be-deficient stars in clus- 
ters younger than ~ 10 Myr could thus provide the genuine 
signature of the accretion process and the proof that proto- 
stars may undergo many extreme bursts of accretion during 
their embedded phases. The search for Be-deficient outliers 
in young clusters, but also in older environments to im- 
prove their statistics, is therefore worth pursuing since their 
existence may improve our understanding of fundamental 
processes in stellar formation and evolution, such as accre- 
tion processes, rotational mixing, and transport of angular 
momentum. 
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